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The two-dimensional surface of a coupled multilayer integer quantum Hall system consists of 
an anisotropic chiral metal. This unusual metal is characterized by ballistic motion transverse 
and diffusive motion parallel (z) to the magnetic field. Employing a network model, we calculate 
numerically the phase coherent two-terminal z-axis conductance and its mesoscopic fluctuations. 
Quasi- Id localization effects are evident in the limit of many layers. We consider the role of inelastic 
de-phasing effects in modifying the transport of the chiral surface sheath, discussing their importance 
in the recent experiments of Druist et al.a 



I. INTRODUCTION 



In a two-dimensional incompressible quantized Hall 
state, the low energy excitations are confined to the edge 
of the sample. These edge states provide a simple way 
to understand transport in both integer and fractional 
quantum Hall systems.B For the integer quantum Hall 
effect with one filled Landau level, there is a single edge 
mode, dcscribablc in terms of a free chiral Fermion. Edge 
states in the FQHE are believed to be (chiral) Luttingcr 
liquids, and have been probedryia tunneling spectroscopy 
in several recent experimentsJj 

In recent years there has been much interest in multi- 
layer quantum Hall systems. In double layer systems the 
layer index plays the role of a pseudo-spin, and these 
systems have revealed a number of new surprises. In 
the opposite extreme with many layers, the samples be- 
come three-dimensional, and a number of new features 
are expected. In such bulk samples with interlayer tun- 
neling smaller than the Landau level spacing, the (in- 
teger) quantized Hall effect in each layer should survive, 
and the sample exhibit a 3d quantum Hall phase. Chalker 
and DohmenQ have recently discussed the phase diagram 
in such a system, in a model of non-interacting electrons 
with disorder. In the absence of disorder, the Landau 
levels will be broadened into bands in the presence of 
interlayer tunneling, t. Disorder further broadens these 
bands. Near the band centers a diffusing 3d metallic 
state is expected. In the tails of the Landau bands, the 
bulk states are localized, but current carrying edge states 
nevertheless lead to a quantum Hall effect. For one full 
Landau level, each layer has a single chiral free Fermion 
edge state, which together comprise a 2d sub-system - a 
chiral surface sheathfl-El This surface phase forms a novel 
2d chiral metal system, which hasJieen analyzed theo- 
retically by a number of authors ElA In the presence of 
impurity scattering, the transport is predicted to be very 



anisotropic, with ballistic in-plane motion and diffusive 
motion parallel to the magnetic field. Vertical transport 
in such a multilayer sample was first investigated exper- 
imentally iruRcf. 12, and has recently been revisited by 
Druist et edn The latter experiment provides striking ev- 
idence of the novel behavior characteristic of the chiral 
metal. 




FIG. 1. Geometry of a multilayer quantum Hall system 
with an interlayer hopping amplitude t allowing vertical 
transport. With N layers the system has height L = Na, 
and a circumference C = 2nr. 

Most of the theoretical work on the chiral metal phase 
has focused on the mesoscopic regime, with the sample 
assumed smaller than the phase breaking lengths. The 
predicted behavior for such phase coherent transport is 
very rich, with three different possible, cegimes (see Fig. |^) 
connected by universal crossovers. Uu'U 

The purpose of this paper is two- fold. First, we revisit 
the phase coherent regime and study in detail the conduc- 
tance and its fluctuations. Performing numerical transfer 
matrix calculations on a directed network model for the 
chiral metal enables us to extract the conductance and its 
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fluctuations in the various regimes, and compare directly 
with earlier analytic approaches. We then address the 
important role of phase breaking processes, which have 
been ignored in earlier theoretical discussions. 

The paper is organized as follows. In Section II we 
briefly review the existing theoretical predictions for the 
phase coherent transport. In Section III we describe the 
network model, and extract numerically the phase co- 
herent conductance in the various regimes. Section IV 
is devoted to a discussion of de-phasing processes, and 
Section V to prospects and conclusions. 



II. PHASE COHERENT REGIME 

For one full Landau level, there is a single free chiral 
Fermion edge mode in each layer, as depicted in Fig. |l|. In 
the presence of an interlayer tunneling amplitude, t (as- 
sumed much smaller than tko c ), these chiral edge modes 
disperse along the z-axis, and form one-half of an open 
2d Fermi surface. Impurities causes electrons to scatter 
about the Fermi surface, as in any dirty metal. Due to 
the chiral nature, the in-plane motion remains ballistic 
with velocity v, even in the presence of impurities. How- 
ever, the (inter-layer) motion parallel to the field becomes 
diffusive, with diffusion constant D. Easier to measure 
than the ballistic velocity or diffusion constant is the z- 
axis (2d sheet) conductivity, a zz , related to v and D via 
an Einstein relations 

D e 2 

o zz = e 2 pD = — • — (2.1) 
va h 

where a is the inter-layer (lattice) spacing and p = 1/ hva 
the density of states. It will be convenient to introduce 
a dimensionless z-axis conductivity via a zz = (e 2 /h)a. 

For a mesoscopic sample with finite circumference, C, 
and number of layers, N = L/a, there are several impor- 
tant time scales. For ballistic motion with velocity v, an 
electron circumnavigates the sample in a time r c = C/v. 
In a time tx = L~/D an electron will diffuse from the 
bottom to the top of the sample. The transport will 
be phase coherent provided the de-phasing time, t^, is 
much longer than both t c and tt,. In principle, this 
mesoscopic regime exists for any sample at sufficiently 
low temperatures, since the de-phasing time diverges as 
T — > (t^ ~ Ti/kftT in the quasi-ld limit of interest). 
Here we focus on the fully coherent regime, returning to 
de-phasing effects in Section V. 

For a sample with finite circumference, C, there are 
two important length scales along the z-axis, which 
demarcate-itue boundaries between three regimes (see 
Fig. ||).EI'B'lI Upon circumnavigating the sample once, an 
electron will diffuse along the z-axis a distance Lq = 
\/Dt c , which can be expressed in terms of the measurable 
z-axis conductivity, a, as, 

L = {acrC) 1/2 . (2.2) 



For finite C with L — > oo the system is one-dimensional, 
and localization along the z-axis is expected. The (typ- 
ical) localization length, £, for such a quasi-ld system 
is proportional to the (dimensionless) Id conductivity, 
£ ~ &idi which can be written, 

£ = 2aC. (2.3) 

Thus both Lq and £ depend only on geometrical param- 
eters, and the measurable z-axis conductivity, a. Notice 
that (£/a) = 2(Lq/o) 2 , so that provided Lo 3> a one has 
t > L . 

As the height L of the sample varies, three regimes 
are possible (see Fig. |J). For L < L <C £, an electron 
typically diffuses from the bottom to the top of the sam- 
ple before circumnavigating the sample once. In this 2d 
chiral metal regime, an electron suffers de-phasing in the 
leads before circling the sample. For L « L « the 
electron circles the sample many times, and phase co- 
herent processes around the sample are important. The 
system behaves like a phase coherent quasi-ld metal. Fi- 
nally, for L 3> £ Id localization effects dominate, and the 
system is a Id (localized) insulator. 




FIG. 2. Three different regimes of the phase coherent 
transport of the surface sheath of a sample with height 
L and circumference C. Both Lo, the typical distance an 
electron diffuses along the z-axis upon circling the sample 
once, and £, the Id localization length, can be deduced 
from the measurable z-axis conductivity and geometrical 
parameters, as discussed in the text. 

The predicted behavior for the phase coherent z-axis 
conductance and it's mesoscopic fluctuations depends 
sensitively on which regime the system is in. Consider 
first the (dimensionless) mean two-terminal conductance 
along the z-axis, G, where the overbar denotes an av- 
erage over disorder realizations. In both the 2d chiral 
metal and the Id metal regimes, ohmic behavior is pre- 
dicted witha, 

G = jo + 0(L/$. (2.4) 

The usual "weak localization" corrections which are of 
order {L/£)° are absent due to the breaking of time re- 
versal invariance. In the Id insulating regime strong lo- 
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calization is operative, and the mean conductance is pre- 
dicted te-fall off exponentially with a universal form (for 

G = 2«/2L)2 exp(-L/2£). (2.5) 

Conductance fluctuations are also of interest, which 
can be characterized by the variance, 8G 2 , where SG = 
G — G. In the|-2d chiral metal and Id metal regimes, 
Gruzberg et aiJiil have shown that the variance can be 
written, 

6G s = $(L/L )+O(L/Z), (2.6) 

where <&(X) is a universal scaling function which 
smoothly connects the two regimes. Deep within the Id 
metal regime the variance approaches a universal number 
well known for quasi-ld metals: Q(L/Lq — > oo) = 1/15. 
In the 2d chiral metal regime, $(L/L ) ~ (Lq/L) 2 for 
L/Lq small. The conductance fluctuations are large in 
this limit, since the sample effectively breaks up into 
t c /t£ = (Lq/L) 2 incoherent regions which add indepen- 
dentl#-jto the conductance and it's fluctuations. Gruzberg 
et alliil have obtained the full universal scaling function, 
§(X), which interpolates between these two limits. In 
the Id localized regime, the conductance is expected to 
be very broadly distributed, with an approximate log- 
normal distribution. 



is the total number of inter-layer tunneling nodes con- 
necting adjacent edge modes (see Fig. ||). The network 
consists of N edge modes, with spacing a and a total 
"height" of L = Na. 

The conductance along the z— axis is obtained by com- 
puting the transmission of electrons from the bottom 
to the top of the sample. Specifically, we use the two- 
terminal Landauer formula to relate the (dirpensionless) 
conductance G to the transmission matrix to 

G = tr[t+t]. (3.2) 



The matrix elements, ty, are the amplitudes for an elec- 
tron incident into channel (or node) i on the bottom edge 
to be transmitted into channel j on the top edge. Here 
N c is the number of channels. 

The transmission matrix is computed numerically by 
iterating a transfer matrix from the bottom to the top 
of the sample. This involves re-expressing each node in 
a form relating the amplitudes in one edge mode to the 
amplitudes in the adjacent edge mode: 

( w ia ) = (ro/to l/to \( v in \ 
V w out J \ l/to r /t J y v out J 



III. NUMERICS 

A. Network model 

Following Chalker and Dohmen,@ we employ a simple 
network model to study phase coherent transport of the 
surface sheath. The network model consists of directed 
links carrying electron current, connected via node pa- 
rameters, as depicted in Fig. ^. All links carry current 
in the x— direction, as appropriate for the chiral surface 
sheath. Scattering at the nodes is characterized by a (real 
and dimensionless) transmission amplitude, to, for tun- 
neling in the z— direction between edge states in neigh- 
boring layers. For a given node the S— matrix relating 
incoming to outgoing amplitudes is given explicitly by, 

( W ° ut ) = (? to )( Win ) , (3.1) 

with t 2 , + r 2 = 1. By construction, this matrix conserves 
the current, |w in | 2 + |w in | 2 = |uw| 2 + |?w| 2 - To model 
the disorder, the electrons are assumed to acquire a ran- 
dom phase along each link connecting adjacent nodes, 
taken for simplicity to be independent and uniformly dis- 
tributed on the interval [0, 2tt]. 

Periodic boundary conditions are taken in the ballistic 
x— direction, with the circumference C = 2bN c , where b 
is the length of a single link in the x— direction and N c 



We study a range of system sizes with the channel 
number N c = 4, 8, 16, 32 and the layer number N = 
8, 10, 12, 16. Being interested in conductance fluctua- 
tions, it is necessary to evaluate the conductance exactly 
for each given disorder realization. The self averaging 
Lyapunov exponents for a sample with L — * oo cannot 
be used to extract the sample to sample fluctuations in 
a finite system. This restriction imposes rather serious 
constraints on the accessible system sizes. 
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FIG. 3. Network model for the surface sheath. Here 
N — L/a — 4 chiral edge modes are interconnected with 
dimensionless tunneling to, with periodic boundary con- 
ditions taken in the x direction of circumference C. The 
z— axis conductance is computed by employing a transfer 
matrix acting in the z— direction. 
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FIG. 4. The single node inside the circle in Fig. 3 
(re-drawn in a more conventional way) represented by the 
matrix in Eqn. (3.3). The transfer matrix progresses from 
bottom to top. 

Since the microscopic parameters of the network 
model, to and b, are not experimentally meaningful quan- 
tities, it is useful to relate them to a macroscopic ob- 
servable, namely the measurable z— axis sheet conduc- 
tivity of the -surface sheath, a. As shown by Chalker 
and Dohmenja this is possible for the network model, by 
summing the Feynman paths analytically. Specifically, 
consider paths that connect the incident electrons on the 
bottom edge to the transmitted electrons on the top edge. 
For C — > oo these paths do not fully circumnavigate the 
sample so that the interference between paths wrapping 
around the sample a different number of times - possible 
for finite C - is completely absent. In the absence of such 
interference the ensemble averaged conductance reduces 
to a sum of classical probabilities: Any two paths which 
pass through a different sequence of directed links will 
have a random relative phase, so that the interference 
term will vanish upon ensemble averaging. To sum the 
classical probabilities of these non- winding paths we fol- 
low Chalker and Dohmenja and consider the transmission 
probability for an electron incident in one channel (say 
i) to be transmitted through N layers: Tn = Ylj l^i| 2 - 
One can then express Tjv+i in terms of Tjv and the sin- 
gle layer transmission probability, Tj 
sum: 



*q, as a geometric 



(3.4) 



where Rn = 1 — Tjv is the reflection probability off N lay- 
ers. Carrying out this geometric sum, gives the recursion 
relation, 



'-N+1 



Ti 



1 

TV 



(3.5) 



which can be readily solved for Tjy. The average con- 
ductance in the absence of interference between paths 
winding around the sample is simply Go = N c Tn, with 
N c the number of channels. It is given exactly by, 



Go 



1 n 



N 1 — tg(l - 1/7V)' 



(3.6) 



as obtained by Chalker and DohmenP for 1/N — > oo. 
The z— axis sheet conductivity follow from Ohms law, 
a = LGq/C, which in the limit G, L — * oo becomes, 



' n 



26 1 - ti 



(3.7) 



Having related the conductivity to the network param- 
eters, the mesoscopic crossover lengths L and £ for a 
finite size netw ork model can be readily obtained from 
Eqn. (|J) and (Oh. 



The exact result for the conductance Eqn. ( |3.6| ) in the 
absence of interference between winding paths should be 
valid even for finite circumference, provided G is large 
enough so that winding paths are rare. The condition for 
the validity of ignoring the interference between winding 
paths is that L <C crC ~ £, so that the sample is in the 
2d chiral metal or Id metal regimes. 

Notice that Go in Eqn. fl3.6| ) is well defined even as 
to — > 1. In this limit, the motion along the z-axis also 
becomes ballistic (for finite TV), and each channel is per- 
fectly transmitted with Go — + N c . It will be convenient 
to define an Ohmic conductance, 



Gohm = Ca/L = £/2L, 



(3.8) 



which coincides with Go when L is large enough that the 
z— axis motion is diffusive. As defined, G h m diverges 
with a as t —* I. The crossover from diffusive to ballis- 
tic z— axis motion occurs when G hm ~ N c . 

The 2d chiral metal regime requires that L -C Lq, or 
equivalently TV <C G h m - However, to avoid a crossover 
into the ballistic regime of the network model requires 
G hm < N c . Thus 2d chiral metal behavior is expected 
for N <C N c . Since this limit is difficult to access numer- 
ically, we focus below primarily on the Id metallic and 
localized regimes. 



B. Results 

In Fig. H we show results for the ensemble-averaged 
two-terminal conductance, G, computed numerically 
from the network model, plotted versus the tunneling pa- 
rameter tp for various different channel numbers, iV c , at 
fixed height, N = 12. The solid lines are the "classical" 



conductance, Go Eqn. (3^), valid in the absence of inter- 
ference between winding paths, and the dashed lines the 
"Ohmic conductance", G h m = Ca/L. Notice that Go 
gives a good fit to the numerical data, except in the low 
conductance regime, G < 1, where Id localization effects 
are expected. The deviations from the classical behavior 
in this regime can be seen more easily in Fig. [|, where 
we plot the same data for the conductance, but now nor- 
malized by Go- Strong deviations are seen for small tg, 
where the system enters into the Id localized regime and 
interference between winding paths is critical. 
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FIG. 5. The mean z— axis conductance G with fixed 
height N = 12 for several different circumferences, plotted 
versus the dimensionless interlayer tunneling probability 
to. The solid lines are Go given by Eqn. (3.6), and the 
dashed lines are the Ohmic conductances Gohm = £/2L 
given in Eqn.'s (3.7) and (3.8). 
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FIG. 6. The mean conductance G from Fig. 5 replot- 
ted after normalizing by Go- The ratio G/Go deviates 
from one as the system becomes localized. 

In order to study the crossover from the Id metallic to 
localized regime, we plot in Fig. [[] the mean conductance 
for N — 12, normalized by G h m — 2£/L, versus 2L/£. 
The data shows a crossover from a Id metallic regime 
with Ohmic behavior, G ~ G hm, to a Id localized regime 
where the conductance vanishes exponentially for L £. 



The solid line is the prediction from Mirlin et al.pQ for the 
mean conductance of a quasi-ld metallic wire obtained 
using supersymmetry methods. The agreement is reason- 
able, but our numerics deviate from the universal form 
of Mirlin et aiJi3 both at large and small £/£. The devi- 
ations at large L/£ are presumably due to lattice cutoff 
effects, since in this regime the localization length along 
the z— axis is comparable to the network model lattice 
spacing a. The deviations for small L/t; are probably 
due to finite size effects. Indeed, as the channel num- 
ber N c increases, the agreement improves. Notice that 
G/G hm vanishes as L/^ — -> (rather than approaching 
unity) due to ballistic behavior in the network model: In 
this limit to — > 1 and G h m diverges whereas G saturates 
at the (finite) channel number N c . 
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FIG. 7. The mean conductance normalized by 
Gohm = £/2L versus 2L/£. For each sample size the 
points correspond to different values of the hopping prob- 
ability to- The solid line is the mean conductance com- 
puted analytically for a quasi-ld system taken from Mirlin 
et al.'s paper. 

In addition to the mean conductance, we have com- 
puted the sample-to-sample conductance fluctuations. 
In Fig. H we have plotted SG 2 versus 2L/£, for height 
N = 16 and various different channel numbers. The 
solid curve is the universal prediction for the variance 
of the conductance of a quasi-ld wire, obtained by Mir- 
lin et alJi3 This curve shows the crossover from the Id 
metallic regime at small where the variance ap- 

proaches the well known universal value, SG 2 = 1/15, 
to the Id localized regime where the fluctuations vanish 
exponentially for L 3> £■ The agreepent between our 
numerical data and the Mirlin et ai.EJ theory is quite 
striking. Again, the deviations for L/£ — » are due to 
the ballistic regime in the network model for to — > 1 
(with finite N) , where the conductance fluctuations van- 
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ish. For L » ^ the localization length approaches the 
lattice spacing. The numerics and theory agree very well 
near the peak in the crossover regime. 
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FIG. 8. Variance of the conductance for different sam- 
ple sizes and hopping amplitudes t$, all plotted versus 
2L/£. The solid line is the variance of the conductance 
for a quasi-ld sample, computed analytically by Mirlin 
et al.. 

Finally, we mention briefly our effort to extract numer- 
ically the conductance in the 2d chiral metal regime. This 
regime requires that L <C L , or equivalently N <§C G hm- 
However, to avoid the ballistic regime when £o — ^ lj we 
must require that G h m < N c , so that we need N N c . 
We have focussed on the conductance fluctuations in this 
regime, since these are predicted to behave very differ- 
ently than in the Id metal, diverging with L/Lq — ► as 
5CP~(L /L) 2 . In Fig. | the variance of the conductance 
is shown for "short" and "wide" samples, with height 
N = 8 and width N c = 16,32,64, plotted versus L/L 
where Lq = yaoC. For each width, N c , we have varied 
the tunneling probability, t§, to get the set of data points. 
The solid line is the analytic prediction from Gruzberg 
et aLjuJ for the conductance variance in the crossover 
regime between the Id and 2d chiral metal. Unfortu- 
nately, the agreement with the analytic result is quite 
poor, although the agreement improves for the widest 
sample with N c = 64. Indeed, the large enhancement in 
the variance for the sample with N c — 64 in the range 
1 < L/Lq < 3 is consistent with the theoretical expec- 
tations. The sharp drop in the conductance fluctuations 
for smaller L/Lq is due to the crossover from diffusive to 
ballistic motion in the network model. The local maxima 
for iV c — 16 at L/Lq « 4 is the same maxima as in Fig. ||, 
and indicates a crossover into the Id localized regime for 
larger L/Lq, where the Gruzberg et al. results do not 
apply. 
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FIG. 9. Variance of the conductance versus L/Lq for 
three different "short" and "wide" samples at various val- 
ues of to- The solid line is the conductance variance com- 
puted analytically by Gruzberg et al. in the universal 
crossover regime between the Id and 2d chiral metals. 
The dashed line is at SG 2 = 1/15 - the value in the Id 
metal regime. 



IV. INELASTIC EFFECTS 

The above results for the phase coherent transport are 
dramatically modified in the presence of phase break- 
ing effects. De-phasing effects can be characterized by 
a phase breaking time, denoted r^, which is the time 
an electron can propagate before having its phase ran- 
domized by interactions with other electrons or phonons. 
In the extreme anisotropic limit of the surface sheath 
with vanishing interlayer tunneling, t = 0, an electron 
propagating in one edge state will interact via Coulomb 
forces with electrons in neighboring edges states, and 
can suffer phase breaking inelastic scattering events. Be- 
ing in Id, the scattering rate, evaluated to leading order 
in the interactions strength u, is linear in temperature: 
1/t0 = c(ulB/2Trhv) 2 kBT/h, with c an order one con- 
stant, Ib the magnetic length and v the edge velocity. 
In practice, the dimensionless ratio uIb/Tt-v is itself also 
of order one, so that 1/t^ ~ ksT/h. For non-zero but 
small interlayer tunneling, the de-phasing rate will prob- 
ably crossover to a two-dimensional T 2 dependence at 
very low temperatures. 

Associated with the de-phasing time are two de- 
phasing lengths: (i) = vt^, the distance an electron 
propagates in the ballistic x— direction before de-phasing 
and (ii) L$ — <J 'Dt^ — y/ aal^ the distance an elec- 
tron diffuses parallel to the field in time r^. Consider 
the transport geometry in Fig. [l|, in which metallic con- 
tacts are applied at z = and z = L. For L^ 3> L, an 
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electron diffuses between the two contacts before being 
de-phased. In this case, transport is mesoscopic, and the 
above phase-coherent results apply. 

For <C L, however, inelastic scattering occurs 
within the sample, and we must reconsider transport 
properties. There are two such important incoherent 
regimes, depending upon the relative magnitude of 1$ and 
G. For 1$ <C C, the electron does not fully circumnav- 
igate the sample before suffering a phase breaking colli- 
sion. In this situation, electron paths which wind a differ- 
ent number of times around the sample do not interfere. 
As a result the system cannot explore the three phase co- 
herent regimes discussed in Sections II and III. Instead, 
the system is appropriately described as a phase incoher- 
ent 2d chiral metal. Nevertheless, there are (small) meso- 
scopic fluctuations expected even in this limit, which we 
discuss below. In the opposite extreme of 1$ ^> C, the 
electron can propagate many times around the sample 
before phase breaking. In this case, the one-dimensional 
motion parallel to the field is phase coherent up to a 
length scale L^. The system should behave like an inco- 
herent quasi-ld wire, with the appropriate (Id) de- 
phasing length, as we discuss further below. 

To describe the transport behavior in these incoherent 
regimes, we employ arguments first applied in Ref. [l^. 
The important observation is that the sample can be sub- 
divided into "patches" , whose size is the maximum area 
over which an electron diffuses in time t^. Each such re- 
gion effectively acts as a classical resistor, and the whole 
sample then as a random resistor network, the properties 
of which are well understood. 

First consider 1$ -C G. Then the patches have di- 
mensions 1$ by Ltj,, and form an array of size C/l$ by 
Denoting by gi the (dimensionless) conductance 
(along the z— axis) of the i th patch, Ohm's law gives an 
average patch conductance of g i = go = al^/L^. The 
conductance fluctuations in each patch, Sgi = gi — go, 
are of order one, - being equivalent to the conductance 
fluctuations of a fully coherent network at the boundary 
between the Id and 2d metal regimes. Since the mean 
conductance can be written go = L^/a, provided the 
patch size is larger than the lattice spacing, L$ ^> a, 
the conductance fluctuations in each patch are much 
smaller than the mean conductance: dgi <C go- In this 
limit, both the total conductance, G, and it's variance, 
SG 2 = G 2 — (G) 2 , can be easily evaluated. A simple 
estimate is to imagine connecting the resistors (patches) 
only vertically (an approximation which gives the correct 
result for the conductance fluctuations up to an order 
one prefactor). Then for each column, the patch resis- 
tances add, so that SG 2 ol « (L^/L) 3 , which is indepen- 
dent of go. Contributing in parallel, the conductances of 
the N co i — C/ltj, columns add, so that the variance of the 
total conductance is simply SG 2 = N co iSG 
be written in the form: 



loi- This can 



SG 2 
G 2 



a 

CL 



a a 
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with G = Ca/L. Notice that the conductance fluctu- 
ations have an appreciable temperature dependence en- 
tering through I growing in magnitude at low temper- 
atures. The mean conductance, however, remains tem- 
perature independent. 

Consider next the Id incoherent limit with 1$ G, 
in which the electron propagates many times around the 
sample before de-phasing. In this limit, the L/ classi- 
cal patch resistors form a one-dimensional random chain, 
and have dimensions G by L$. Due to Id localization ef- 
fects, the conductance of each such segment will depend 
strongly on it's length, L$, and hence on the tempera- 
ture T . For example, when is much smaller than the 
Id localization length £, the (mean) conductance of each 
segment is given by, 

G seg (L ) = (aC/L^) - + O(L^/0 2 , (4.2) 

where the first term is Ohm's law, and the second term 
reflects the leading Id localization corrections within the 
unitary ensemble. In the opposite limit, 3> £, one ex- 
pects a stronger length (and temperature) dependence, 
G S eg(-k</>) ~ ex p( — Ls/2£). The total conductance fol- 
lows by simply adding the series resistances of each of the 
L/Lcf, segments. In the Id metallic regime with L$ <§; £, 
this gives, 



-=j _ aC 2 L$ Lff, 



(4.3) 



which depends on temperature through L^(T). 

Experimentally, such conductance fluctuations are 
usually observed not by looking at different samples, but 
by varying the applied magnetic field in such a way as to 
change the phases accumulated by interfering electrons 
and thereby effectively change the disorder. The conduc- 
tance fluctuations in this context are characterized not 
only by their amplitude, discussed above, but also by a 
characteristic field scale SB^. This scale is defined by the 
amount the applied field must be changed in order that 
the conductance of a fixed sample becomes uncorrelated 
with its previous value. Physically, the conductance fluc- 
tuations arise from constructive interference of two paths 
enclosing an area of the phase-coherent patch size. The 
total change in phase shift around this loop in units of 
2tt is simply the change in magnetic flux through this 
area divided by the flux quantum (f>o — hc/e. The char- 
acteristic field 8B,p, which changes the phase around the 
loop by 0(tt), is thus simply the field which puts, say, 
half a flux quantum through this coherent area. Assum- 
ing the magnetic field has a non-negligible angle to the 
surface sheath (which we believe to be the case in the 
experiments of Druist et ah), this gives 



SB d 



4>o/l<pL<p l<f 
4>ojCL^ l 4 



< G 
> G 



(4.4) 



(4.1) 



in the two incoherent regimes. Note that since 1$ and 
increase as temperature is lowered, the conductance 
varies very rapidly with field at low temperatures. 
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V. CONCLUSIONS 

We conclude with a comparison of these theoretical re- 
sults to the experimental data of Druist et alE. Druist et 
al. have measured the z-axis transport in a series of mul- 
tilayer quantum Hall samples. Specifically, the samples 
consisted of 50 layers of 150A GaAs layers alternating 
with 150A Alo.iGao.gAs barriers doped at their centers 
with Silicon. The vertical separation between each of 
the 50 2d electron gases is a — 300A A simple Kronig- 
Penney analysis gives an estimate for the z-axis band- 
width of t = 0A2meV. When the applied magnetic field 
was tuned onto an integer quantum Hall plateau, the 
z-axis conductance - dropping with temperature - was 
found to saturate below about 200m_ft'. Since the low 
temperature z-axis conductance scaled linearly with the 
circumference (perimeter) of the samples, which ranged 
from 400/im < C < 7mm, Druist et. al. argued that the 
conduction was being dominated by the 2d chiral surface 
sheath. The resulting sheet conductivity along the z-axis 
was found to be a « 4 x 10~ 4 on the v — 1 plateau, and 
about a factor of three larger for v = 2. 

A theoretical estimate for the z-axis conductivity of the 
surface sheath at one full Landau level can be obtained 
fromO 



al g 



h 2 v 2 



(5.1) 



where Iq is an elastic mean free path for edge scattering 
and v is the (ballistic) edge velocity. Unfortunately, both 
v and Iq are difficult to estimate reliably, depending on 
the detailed slope and irregularities of the edge confining 
potential. However, we expect that in the limit of large 
magnetic field, Iq > Ib, where Ib is the magnetic length 
(lo may grow much longer than Is as the edge is made 
cleaner). Moreover, we expect v to be bounded above 
by the edge velocity for a hard-wall confining potential, 
so that v < w c /s/27r, with uj c the cyclotron frequency. 
Putting in these (rough) bounds, we obtain 



> 



(27r) 2 t 2 a 
h lu?.Ib 



(5.2) 



Using the parameters appropriate for the Druist et al. 
experiment, this gives a > 6 x 10~ 5 , about an order of 
magnitude smaller than the experimental value. Given 
the uncertainties in v and Iq, as well as possible shifts 
in t due to interaction effects, this level of agreement is 
reasonable. 

Taking now the measured value of a, we can estimate 
the two length scales which determine the system be- 
havior in the mesoscopic limit. The samples studied by 
Druist et al. had a range of circumferences 400/im < 
C < 7mm, which correspond to length s 2 < Lp/ g < 10 
and 10 < £/a < 200, upon using Eqns. dU)-(|J). Since 
N = L/a = 50 in these experiments, in the mesoscopic 
limit these samples should span the quasi-ld metal and 



Id localized regimes. At low temperatures, we would 
therefore expect a strong suppression of the conductiv- 
ity and significant temperature and circumference depen- 
dence, especially in the smaller samples. That such ef- 
fects are not observed must be attributed to inelastic 
effects. Indeed, as shown below, estimates for the in- 
plane de-phasing length 1^ give 1$ <C C even at the lowest 
temperatures and for the smallest sample. In this limit, 
mesoscopic effects are greatly suppressed, and the system 
is best thought of as an incoherent 2d chiral metal. This 
accounts naturally for the observed low temperature sat- 
uration of the conductivity (it remains to be seen whether 
the weak residual temperature dependence at low T can 
be fitted to the expectedB form cr(T) — c(0) cx T 2 ). 

We can attempt to estimate the de-phasing length 1$ 
via 



A 



h u 
uIr 



hv 



(5.3) 



however there is considerable uncertainty in the param- 
eters - particularly the edge velocity v. As a crude esti- 
mate we take A = 1, a dimcnsionless interaction strength 
of unity uls/hv = 1 and an edge velocity estimated for 
a hard- wall confining potential v — uj c Ib/2t:. In the 10 
Tesla field used by Druist et al. in the v — 1 plateau and 
at the lowest temperatures studied of T — 50mK this 
gives the rough estimate Ij, ~ 20/im 

Fortunately, one can also extract estimates for 1$ di- 
rectly from the experimentally measured conductance 
fluctuations. In fact, this can be done in two ways, 
thereby providing a consistency check. One determina- 
tion i s fro m the amplitude of the fluctuations. Solving 
Eqn. (LI) gives 



G 



-(5G Z 



(5.4) 



Because the fourth power of 5G appears above and the 
amplitude A is unknown, there is again considerable un- 
certainty in l^. For the Druist et al. experiments, we 
obtain l^ k, 26/im, consistent with the above theoretical 
estimate. 

A second determination comes from the magnetic field 
scale of the conductance fluctuations. From the above 
estimates, we see that Ltf, = yjcral^ < a (using the mea- 
sured a = 4 x 10~ 4 ). This is close to the "incoherent 
tunneling" limit, and we e xpect it is appropriate to re- 
place — > a in Eqn. (4.4), giving 



(5.5) 



For the Druist et. al. experiment, this gives l^ « 3/xm at 
T = lOOmAT, somewhat smaller than the first estimate. 
In this case there are also considerable uncertainties due 
primarily to incomplete knowledge of the degree of inter- 
layer flux penetration. However, all three of the above 
estimates give 1$ <C C. 
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In summary, the experiments so far are consistent with 
the picture of an incoherent 2d chiral metal. Several op- 
portunities exist for further theoretical and experimen- 
tal study. Samples with smaller circumferences in the 
range of 10 to 20^m would be highly desirable, since the 
mesoscopic regime would then be accessible below several 
hundred mK. In this limit, the rich and varied crossovers 
between the three mesoscopic regimes could be accessed 
experimentally. Theoretically, a more quantitative study 
of inelastic scattering and de-phasing lengths would be 
desirable in order to achieve a precise comparison with 
experiment. Particularly interesting from both points of 
view is the temperature dependence of which we be- 
lieve should exhibit linear scaling with temperature over 
a broad range. A field-theoretic treatment of de-phasing 
effects could be useful in providing the desired tighter 
link with experiments. 
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